242 


SPECTRAL ANALYSIS OF THE FORBUSH DECREASE 
OF JULY 13, 1982 


SH 5.1-7 


E. Vainikka, J.J. Torstl, E. Valtonen, H. Lumme, 

M. Nieminen, J. Peltonen, and H. Arvela 

Department of Physical Sciences, University of Turku, and 

Wihuri Physical Laboratory, University of Turku 
SF-20500 Turku, Finland 

ABSTRACT 

The maximum entropy method has been applied in the 
spectral analysis of high-energy cosmic-ray intensity 
during the large Forbush event of July 1 3 , 1982 . An 
oscillation with period of about 2 hours and amplitude 
of 1 - 3 Z was found to be present during the decrease 
phase. This oscillation can be related to a similar 
periodicity in the magnetospheric field. However, the 
variation was not observed at all neutron monitor 
stations. 

In the beginning of the recovery phase, the intensity 
oscillated with a period of about 10 hours and ampli- 
tude of <,3 

1. INTRODUCTION 

One of the largest Forbush decreases (FDs) was observed on July 13, 
1982. It was caused by a solar flare of type 3D on July 12 at 09 1 6 UT 
at position 11 N, 37 E. The arrival of the shock front was indicated 
by a sudden commencement (SSC) at 1 6 1 8 UT on July 13. The magnetic dis- 
turbance generated the largest geomagnetic storm ever observed In 
Finland. At Nurmi jSrvi Geomagnetic Observatory (60.5 N, 2*1.7 E) , the 
horizontal component of the geomagnetic field reduced 20 $ within a few 
hours, and returned to its predecrease level in 12 hours. 

The decrease of high-energy cosmic rays began at about 1620 UT on 
July 13, and the recovery phase lasted several days. The amplitude of the 
decrease was of the order of 20 % as measured by the Turku double neutron 
monitor and other high-latitude monitors. 

In order to find out possible periodicities In the cosmic-ray Inten- 
sity during thi s FD, we applied the Maximum Entropy Method (MEM) , proposed 
by Burg In 1967, In analyzing data from nine neutron monitor stations. 

The advantage of MEM Is that It leads to much better resolution than the 
conventional Blackman-Tukey method. It also gives more realistic power 
spectral estimates, especially for short data records. 

2. DATA AND METHOD OF ANALYSIS 

The neutron monitor data used In our analysis are displayed In Figure 
1. Before the frequency analysis, long-term trends (period > 12 hours) 
were removed by using the moving average method. 

The power spectral density, P(f), was estimated from the equation 

P(f) - At P M+1 |l + E a M|< exp(-2rci fkAt) | " 2 , 

k«1 
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— 9‘ . 1 ‘ N< - Utr °n monitor counting rates on July 12-15, 1982 ai 
stations- (see table i for the abbreviations) .The l^est 

iarv?°^f??„! '. co '»P°nent of the geomagnetic field as measured 
jarvi station in Finland. 


nime 
shows 
at Nurmi- 


where f is the frequency, At is the ti 

a " M ,en 9 th °f the prediction 

the coefficients a Mk were determined f 


me interval between the data points, 
error filter. The constant P and 
rom the matrix equation M+1 



where $ ( 1-0 n ) is 

of the method, see e.g. 
5 or 15 min, and M was 
of data points). 


the autocorrelation with time 
Ulrych and Bishop 1975). The 
determined in the standard way 


lag nAt (for details 
time interval .At was 
(M4N/3, H ■ number 


3. RESULTS AND DISCUSSION 

3.1. Decrease phase. The decrease phase started at about 1620 UT. However 

th n ° nS ? t tlme Was about ,i<30 UT - The explanation Is that the ’ 
asymptotic directions of Inuvtk rotate with the Earth so that at 1^*30 UT 
nuv receives particles from the direction of the Interplanetary 
magnetic field-line. But these particles were strongly scattered by the 
arriving shock front, and thus their intensity decreased earlier. 

e , on ^ ft cant variation In the frequency range 3*10*5, IQ-3 Hz 

occurred with periods between 112 and 126 min (Table 1)7 The amplitude 
was largest at the mountain station Jungfraujoch (*3 %, as reported 
also by Debrunner-et al. 1 983) . At Finnish stations (Turku and^ulu) the 
amplitude was £2 %. Aq analysis of the neutron-multiplicity data of the 
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Table 1 . The most significant periods in the cosmic-ray intensity at 
various neutron monitor stations. 

Period (min) 

Station R vert (GV) Decrease phase Recovery phase 


Jungf raujoch 

(Ju) 

4.5 

118 

716 

Kiel 

(Ki) 

2.3 

115 

706 

Turku 

(Tu) 

1.1 

115 

706 

Oulu 

(Ou) 

0.8 

112 

788 

Hermanus 

(He) 

4.6 

- 

430 

Deep River 

(De) 

1.1 

- 

640 

Goose Bay 

(Go) 

0.6 

- 

596 

1 nuvi k 

(In) 

0.2 

126 

540 

Alert 

<A1) 

0.0 

(119) 

702 



Turku neutron monitor showed that this vari- 
ation was no more present at multiplicities 
24 corresponding to median primary rigidities 
*35 GV. 

The power spectral densities are shown in 
Figure 2. It is noteworthy that the 2-hour 
periodicity was practically absent at Hermanus 
Deep River, and Goose Bay. Further, the period 
119 min of Alert is not statistically signifi- 
cant. 

In Figure 3, the power spectral density 
is presented in more details around the 2-hour 
peak. The superiority of MEM to the conven- 
tional Blackman-Tukey method is evident. 

The origin of these oscillations is clear- 
ly in the strong variation of the magneto- 
spheric field as can be seen in Figure 4, 
where the observations of the satellite GOES-2 
are reproduced. 


3.2. Recovery phase . In the beginning of the 
recovery phase, rather strong oscillations 
were observed at all stations. The only re- 
markable variations occurred with a period of 
about 10 hours and amplitude of 3 % (Table 1 
during the decrease phase, an( j pjg Ure 5), a similar periodicity was found 


Fig. 2 . Power spectral 
densities of the cosmic- 
ray intensity observed 
at various stations 





Fig* 4 . The total magnetic field as 
measured by the satellite GOES-2 
(from Parsignault et al . 1983). 

PERIOD (h) 

the shock front. Power spectral densities 

Inspection of Figure 1 reveals durin 9 the recovery phase, 

some phase shifts of a few hours in 

the 10-hour variation (e.g. Deep River and Goose Bay vs* European sta- 
tions). This feature Is understood In terms of differences in the 
asymptotic directions of the stations. 

4. CONCLUSIONS 

The unusual periodicity of 2 hours during the decrease phase of 
the FD of July 13, 1982, Is a direct consequence of strong oscillations 
in the magnetospher i c field. The mechanism of the interactions between 
the solar wind and the Earth's magnetosphere leading to these oscilla- 
tions requires further studies. 
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